Noisy nonlinear dynamics of vesicles in flow 
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We present a model for the dynamics of fluid vesicles in linear flow which consistently includes 
thermal fluctuations and nonlinear coupling between different modes. At the transition between 
tank-treading and tumbling, we predict a trembling motion which is at odds with the known deter- 
ministic motions and for which thermal noise is strongly amplified. In particular, highly asymmetric 
shapes are observed even though the deterministic flow only allows for axisymmetric ones. Our re- 
sults explain quantitatively recent experimental observations [Levant and Steinberg, Phys. Rev. 
Lett. 109, 268103 (2012)]. 

PACS numbers: 87.16.D-, 47.63.-b, 05.40.-a 



Introduction.- For a system close to a dynamical insta- 
bility, the effect of small perturbations can be dramatic 
[l|42(. In particular, external noise [2| such as thermal 
fluctuations may be strongly amplified in the vicinity of 
a bifurcation, therefore crucially affecting the dynamics 
of the system. Such effects were observed, e. g., in binary 
fluids [4|, fiber lasers [5| and even population dynamics 
[fj. Recently, Levant and Steinberg j7( realized exper- 
iments in which such an amplification of thermal noise 
occurred for fluid vesicles subject to an external flow. 

Fluid vesicles [8| are closed lipid bilaycr membranes 
which mimic biological objects such as red blood cells. 
Understanding their dynamics in flow may therefore help 
to elucidate blood rheology [9l— ll 3| or chemical activity 
14| . This dynamics in strongly nonlinear due to curva- 
ture elasticity and the incompressibility of the fluid mem- 
brane. Most theoretical models 115142(1 and numerical 
simulations 
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-24( neglect thermal fluctuations, mainly 
because the bending rigidity of the membrane is of the 
order of k ~ 2bk-oT. For weak flows obeying Stokes equa- 
tion, these models predict three main dynamical regimes: 
(i) a tank-treading motion (TT), for which the orienta- 
tion of the vesicle remains constant while the membrane 
rotates around the interior fluid; (ii) a solid-like rotation 
of the vesicle called tumbling (TB); (iii) and a vacillating- 
breathing (VB) motion 251 126| (also called swinging [27[ 
or trembling [16() at the TT-TB transition, in which the 
orientation angle oscillates around whereas the shape 
periodically deforms. An important feature of these pre- 
dictions is that the vesicle shape remains axisymmetric 
at all times. 

The TT and TB motion observed in experiments |28l — 
|30| are well described by the theoretical predictions. 
However, a series of experiments (71. I29L l3lM33j | shows that 
the trembling motion (TR) found at the TT-TB transi- 
tion is very different from the deterministic VB motion. 
The experimental TR is an aperiodic motion character- 
ized by highly asymmetric shapes and a larger basin of 
attraction in the parameter space than VB. The discrep- 
ancy between model and experiments is presumably due 
to thermal fluctuations [7|, 133J , a claim supported by some 



qualitative results obtained in stochastic simulations [3J 
|36| which, however, have not looked systematically at the 
TT-TB transition. Analytical models including thermal 
fluctuations 37H39J| exist but have not yet addressed the 
TR motion. 

In this Letter, we present a model for vesicles in planar 
linear flow which systematically and consistently includes 
thermal fluctuations. We predict a trembling regime at 
the TT-TB transition which is at odds with deterministic 
predictions but agrees very well with the experimental 
observations. We show that, because of the amplification 
of thermal noise, highly asymmetric shapes are observed 
even though the deterministic external flow would only 
allow for axisymmetric ones. 

Model- We consider a vesicle of volume V = 47t7?q/3 
and surface area A = (4ir + A)7?q, where Rq is the radius 
of a sphere of the same volume and A the excess area 
relative to it. This vesicle encloses a fluid of viscosity rji 
and is suspended in another liquid of viscosity r\ . The 
viscosity contrast A = r\ijr\ is set to 1 in the following. 
The two-dimensional fluid membrane is assumed to be 
incompressible and impermeable so that V and A, or 
equivalently Rq and A, remain constant. This vesicle is 
subject to a planar linear flow, for which the unperturbed 
velocity field can be written as 



v = s(xe y + ye x ) + uj(ye x - xe y ) 



(1) 



where s is the strength of the extensional flow, which 
tends to stretch the vesicle along the (y=x)-direction, 
and u> the vorticity, which is the strength of the rotational 
component of the flow. A pure shear flow corresponds to 
s = u). Wc consider weak flows such that the Reynolds 
number is low and the hydrodynamics is described by 
the Stokes equation. Moreover, the energetic cost of the 
vesicle deformation is given by the Helfrich energy [8| 



H= dA 
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(2) 



where k is the bending rigidity of the membrane, H the 
local mean curvature and a is the surface tension ensur- 
ing local and global incompressibility. 



In order to derive analytically the equations of motion, 
we assume that the vesicle is quasispherical, i. c, its 
excess area A is small. We then expand the radius 
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in the basis of spherical harmonics 3^ m . Adapting |37| , 
we solve the hydrodynamic problem inside and outside 
the vesicle and match the solutions at the sphere of radius 
i?o, which leads to the equations of motion 



ui m = vmuuim - Vi 



an 

du* 



- ihsSi t2 (6 m , 2 - 5 m ,-2) + (lm(t) 

(4) 
for the coefficients ui m , where the dot stands for the time- 
derivative and h = 24-\/27r/15/Tl, /.t; = Ti/t] Rq and 
T ; = 1(1 + l)/(4/ 3 + 6l 2 - 1). The thermal noise Q m (t) is 
white and its correlations 

(Cim(*)Ci'm'(0> = 2mk B T(-l) m 6 l , v 6 m ,- m >6(t-l/) (5) 

are given by the fluctuation-dissipation theorem in equi- 
librium. The first term on the right-hand-side of Eq. 
(|4]) is a norm-preserving rotation of ui m for all harmon- 
ics and is due to the vorticity w of the external flow. 
The third term is only present for the prolate harmon- 
ics 3^2, ±2, which are the only ones directly excited by 
the flow. For this reason, wc can distinguish between 
different regimes by looking at the evolution of the ori- 
entation angle 022 of the prolate harmonic defined by 
U22 = |w22|e _2l<ji22 . It corresponds to the inclination an- 
gle of the long axis of the vesicle in the flow plane. Its 
equation of motion follows from ((4|) as 
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where (^ is a Gaussian white noise having the correla- 
tions ([S]). We see that for small values of the ratio lj/s, 
the deterministic part of Eq. ^ may have a fixed point, 
whereas for larger vorticities, the right-hand side is al- 
ways negative. These two cases correspond roughly to 
the TT and TB motions, respectively. 

The nonlincarity of the equations ([?]) is due to the term 
depending on the Hclfrich energy ©. Following Rcf. 
|16j |. we expand the bending energy T-L to third order in 
the coefficients ui m [40( . Such terms are necessary for the 
dynamics to be insensitive to initial conditions 17|. We 
do not include the hydrodynamic corrections discussed 
in Refs. [18|, |41|, but keep in mind that they become 



relevant for strong enough flows. Through the condition 
of constant volume, we can express «oo as a function 
of all other coefficients with I > 1. The homogeneous 
part of the surface tension a has to be determined such 
that A remains constant. The local incomprcssibility of 
the membrane has already been taken into account in 
deriving Eq. (fj[|). 
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FIG. 1. Temporal evolution of the inclination angle i/>22 (solid 
black curves) for A = 0.3, S = 30 (s = 0.032 s" 1 ) and A = 0.8, 
1.4, 1.5 and 1.8 (u> = 0.066 s" 1 , 0.11s" 1 , w = 0.12 s -1 and 
to = 0.15 s" 1 ). The blue dashed curves are the corresponding 
deterministic trajectories. 



In the deterministic case, only the harmonics of order 
I = 0, 2 and 4 survive to that order of expansion [18J]. In 
particular, terms with odd / remain strictly equal to zero. 
Now, if wc include thermal fluctuations, all ui m have a 
positive mean absolute value. Therefore, they cannot 
be discarded out of symmetry considerations since they 
might be important at the dynamical transitions. We 
will thus keep all harmonics between I = 2 and I = 5 
(the modes with I = 1 correspond to a translation of the 
vesicle as a whole and are not relevant to our problem) 
and solve numerically the eighteen coupled equations of 
motion. In the numeric scheme, we replace the surface 
tension a, which strictly enforces the constraint of con- 
stant area, by a softer energetic constraint of the form 
K(A^ — A) 2 /2 in the bending energy (J2|), where ^( 3 ) is 
the expansion of the area A to third order in Ui m . For 
K — > oo, we thus recover the same dynamics as with the 
rigid constraint a. The soft constraint, however, is much 
easier to treat numerically since it avoids the intricacies of 
multiplicative noise [421 ] . Furthermore, in order to com- 
pare our results to other models and to experiments, we 
introduce the parameters [iff 



s ^ Uir ri Rl s 
~ 3^3 «A S 
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where S is the dimcnsionlcss strength of the clongational 
flow and A is the rescaled ratio of the rotational and 
clongational components of the flow. It has been claimed 
that these two quantities are enough to describe the dy- 
namics of the system [l6( , a fact contested in other the- 
oretical models |43( but found to be pretty accurate in 
experiments [31(. In the following, we set T = 293 K, 
Rq = 15 pm, k = 25&bT and r\ = ImPas, which corre- 
spond to typical experimental values [7| . The comparison 
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FIG. 2. Top: evolution of some mode amplitudes for A = 
0.66, S = 60.7 and A = 1.8 (corresponding to s = 0.14s" 1 
and oj = 0.44s" 1 ). Bottom: instantaneous shape contours in 
the flow plane at four different times (see vertical lines in the 
plot). 



to the deterministic case is done by ignoring the thermal 
noise in our equations. 

Regimes of motion.- Fig. [1] shows the temporal evo- 
lution of the inclination angle ^22 for S = 30 and four 
different values of A. For small values of the vorticity 
(A = 0.8), the vesicle is in the TT regime, in which its 
orientation fluctuates around the deterministic value but 
remains positive. If we increase the vorticity of the flow, 
we observe the TR motion (A = 1.4 and 1.5), for which 
the oscillation angle irregularly oscillates around 0, al- 
though some full rotation might occur from time to time 
(see beginning of A = 1.5). Without thermal noise, the 
model predicts a periodic VB motion for A = 1.4 but 
TB already for A = 1.5, whereas the stochastic motion 
is still TR. This observation is in agreement with exper- 
imental results, which find that TR happens for a much 
broader range of parameters than predicted by determin- 
istic models 3l|, |32|. For A = 1.8, we see a mixture be- 
tween TR and TB, with bursts of TR happening between 
the TB cycles. Thermal fluctuations mix both regimes at 
the TR-TB transition, a similar effect to what happens 
for undcformable capsules at the transition between tum- 
bling and swinging [39(. 

The difference between stochastic and deterministic 
models becomes even clearer if one looks at the evolu- 
tion of the shape of the vesicle. Fig. [2] shows the time- 
evolution of the coefficients M20 and uu (I = 2. ..5) for a 
typical mixed TR/TB motion, using the same parameters 
as in Fig. 2 of Ref. [7J. We choose these coefficients be- 




FIG. 3. Power spectra of the rescaled amplitudes \uu\ 2 /A 
as a function of the mode number I = 2. ..5 for S = 30 and 
A = 0.8 (TT, green dash-dotted), 1.5 (TR, red full) and 2.5 
(TB, black dashed) and four different excess areas A. The 
blue dotted lines are the corresponding equilibrium spectra 

a- 



cause they determine the shape of the vesicle in the flow 
plane, which is observed in experiments. In the deter- 
ministic case, the odd harmonics U33 and M55 are strictly 
equal to zero, whereas 1*44 remains much smaller than 
M20 and M22 ■ The shape of the vesicle remains ellipsoidal, 
even though there is a periodic "breathing" due to the 
interplay between U20 and «22- With thermal noise, all 
modes are excited and higher harmonics become regu- 
larly larger that the second one, leading to strongly dis- 
torted shapes as shown by the four snapshots of the con- 
tour of the vesicle in the flow plane, which are similar to 
the experimental ones [7( . This shape distortion happens 
when the inclination angle becomes negative, i. e., when 
the vesicle is compressed rather than stretched by the 
external flow. This phenomenon is reminiscent of vesi- 
cle wrinkling [4J-|46j, which happens when the direction 
of the elongational flow is suddenly inverted. However, 
wrinkling is a transient, deterministic phenomenon which 
is triggered by a discontinuity in the external flow. The 
trembling motion presented here exists only with thermal 
fluctuations and is a steady-state property of the vesicle 
dynamics. Note that this kind of shape deformations are 
also present in TB, since the vesicle also passes through 
negative inclination angles, but become smaller for larger 
A because the vorticity ui and thus the rotation frequency 
of the vesicle become larger. 

Power spectra and phase diagram.- We can quantify 
the shape dynamics by looking at the power spectra of 
the rescaled harmonics |u/;| 2 /A in the different regimes, 
as shown in Fig. [3] for S = 30 and different values of A. 
For low excess areas (A = 0.1), the power spectra in TR 
and TB are very close to the equilibrium one due to the 
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FIG. 4. Mean velocity < 022 > rescaled by the deterministic 
one < 022 >dot as a function of the dimensionless variables S 
and A for A = 0.1 and 0.3. The white lines correspond to the 
TT-TR transition. The four crosses correspond to Fig. [l] 



fact that the absolute strength s of the elongational flow 
is relatively small (s ~ 0.01s -1 ) thus exerting stresses 
comparable to the thermal ones. The TT spectrum al- 
ready shows a larger second harmonic than in equilib- 
rium (and correspondingly smaller higher-order harmon- 
ics) because the orientation angle 022 remains positive 
and therefore the vesicle is constantly stretched. S being 
constant, larger A means larger s, which should suppress 
the effect of thermal fluctuations. This suppression in- 
deed happens for the TT and TB motions, for which we 
see a decrease of the amplitude of the odd modes while 
the fourth mode, which is coupled to the second one, be- 
comes larger. By contrast, in the TR regime the third 
mode remains always one order of magnitude larger than 
the fourth one even for large A. These features are also 
found in the experimental data [7], l33[ although the cou- 
pling between even modes seems stronger there. The TR 
motion is thus directly caused by the strong amplification 
of thermal noise at the dynamical TT-TB bifurcation. 

Further interesting features are revealed by the dynam- 
ical phase diagram of the system, as shown in Fig. [4] for 
two different excess areas. Since the transitions between 
the different regimes are smoothed out by thermal fluctu- 
ations, we choose the mean tumbling velocity < 022 > as 
an order parameter. When the deterministic prediction 
is TB, we rescale this order parameter by the determin- 
istic mean velocity < 022 >dct> so that we clearly see 
the effect of thermal noise. In the TB regime, this ratio 
is close to 1 (white regions) while it is almost in the TT 
and TR regimes (black regions). The gray regions indi- 
cate a mixed TR/TB motion (see Fig. [Qfor A = 1.8). We 
further define TT as the regime in which the mean ori- 
entation angle is larger than its standard deviation, such 
that we are able to distinguish between TT and TR. For 
small enough S, the dynamics is dominated mainly by 
diffusive motion of the orientation angle, making it hard 
to distinguish between different regimes. For this rea- 
son, these phase diagrams do not extend to S = 0. The 
first striking feature is that the TR region is much larger 
for A = 0.1 than for A = 0.3. This difference is due 
to the fact that, for fixed S, the magnitude of the elon- 



gation s grows with A, such that thermal fluctuations 
have a greater effect for small excess areas. The TR re- 
gion is also much wider than the one predicted by the 
deterministic model of [16J , which corresponds approxi- 
mately to our model without noise. Going to larger A, 
the TR region becomes thinner, but then the neglected 
higher-order hydrodynamic terms become important. If 
we had included these, we might observe a broadening 
of the TR region with growing A, similar to what hap- 
pens for the VB region in deterministic models including 
hydrodynamic corrections [18|. In addition, the TT-TR 
transition happens approximately at fixed values of A for 
all A. The phase diagrams for different values of A have 
therefore a large overlap, which might explain why S and 
A are so effective for classifying experimental data [33J. 
Our results show, however, that A is also important for 
the dynamics and that S and A are not sufficient control 
parameters. Finally, because of the thermal fluctuations, 
the size of the TR region is quantitati vely comparable to 
the one measured in experiments [3ll-l33[. 

Conclusion.- We have shown that even for large vesi- 
cles with an effective radius of 15 um, thermal fluctua- 
tions are crucial to understand the dynamics in linear 
flow. This is due to the sensitivity of the nonlinear dy- 
namics of the vesicle to small perturbations close to the 
TT-TB transition. In this case, the small perturbations 
corresponds to the always present thermal fluctuations 
which incessantly destabilize the system, generating the 
trembling motion. Our results thus give an explanation 
for the strong shape fluctuations observed in experiments 
that are not predicted by deterministic models. The gen- 
eralization of these results to other soft objects such as 
red blood cells or elastic capsules will be the focus of 
further investigation. 

We thank V. Steinberg for discussions and acknowledge 
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